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Abstract 

We study various channels of top and Higgs Flavor Changing Neutral Couplings 
(FCNC) in the two Higgs doublet Model with natural flavor conservation (2HDM-I 
and 2HDM-n). We update known results about t — > 07 , cZ , eg and comment on t — > 
ch^. The decays t — > ch^ as well as {h^, H^} — > tc are sensitive both to the bottom 
Yukawa coupling as well as to the trilinear scalar couplings h?H^H~ and H^H~^H~ . 
After imposing unitarity constraints as well as vacuum stability conditions on scalar 
sector parameters, in 2HDM-n we found that for large tan /3 > 40 and rather light 
charged Higgs mass Mh± < 150 GeV, the maximum values allowed for Br{t — > ch^), 
Bi{H^ ic) and Br(/i° tc) are: 8 x 10~^, 10~^ and 10~^ respectively. For charged 
Higgs mass in the range [200, 300] GeV, which can accommodate B — > Xsj constraint 
if one takes into account large theoretical uncertainties, the branching ratio of both 
tc and — > tc can still be slightly larger than 10~^. For — > tc, its 
branching ratio is smaller than ~ 10~^ in both 2HDM-I and 2HDM-II. We study 
also the top-charm associated production at e+e~ colliders and its 77 option as well 
as at muon colliders. It is found that the cross section of 77 — > tc can be of the 
order 0.01 — > 0.1 fb near threshold region, while the cross section of e^e~ — > tc is 
well below 10~^ fb. The situation is slightly better for muon colliders where a few fb 
cross section can be reached for large tan /? and low center of mass energy < 500 
GeV. 
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One of the goals of the next generation of high energy colhders, such as the large hadron 
colliders LHC [1] or the International Linear Colliders ILC [2] or muon colliders [3], is 
to probe top flavor-changing neutral couplings 'top FCNC as well as the Higgs flavor- 
changing neutral couplings 'Higgs FCNC. FCNC of heavy quarks have been intensively 
studied both from the theoretical and experimental point of view. Such processes are being 
well established in Standard Model (SM) and are excellent probes for the presence of new 
physics effects such as Supersymmetry, extended Higgs sector and extra fermions families. 

Within the SM, with one Higgs doublet, the FCNC Ztc vanish at tree-level by the 
GIM mechanism, and the ^tc and gic couplings are zero as a consequence of the unbroken 
SU{?>)c X U{1)q-^ gauge symmetry. The Higgs FCNC Htc and Hhs couphngs also vanish 
due to the existence of only one Higgs doublet. Both top FCNC and Higgs FCNC are 
generated at one loop level by charged current exchange, but they are very suppressed by 
the GIM mechanism. The calculation of the branching ratios for top decays yields the SM 
predictions [4], [5]: 

Br(t ^ Zc) = 1.3 X 10"^^ Br(t ^ 7c) = 4.3 x 10^^^ Br(t ^ gc) = 3.8 x 10-^\ (1) 
Br(t ^ He) = 5.6 3.2 x 10~^^ for Mh = 115 ^ 130GeV. (2) 
Br(H ^ tc) fti 1.5 X 10"^^ (resp 3 x IQ-^^) Mh = 200 (resp 500) GeV (3) 

In the last decade, there has been intensive activities to explore the top FCNC both 
from the experimental and theoretical point of view. Experimentally, CDF, DO [6] and 
LEPII [7] collaborations have reported interesting bounds on top FCNC. These bounds 
are very weak but will improve in the next years, first with Tevatron Run II, and later 
with the future experiments such as LHC and/or ILC. 

From the theoretical side, many SM extensions predict that these top and Higgs FCNC 
can be orders of magnitude larger than their SM values (see [8] for an overview). Among 

those extensions, we cite: 

i) Models with exotic quarks [9], [10] where one can have branching ratios of the order 
10"^, 10~^ and 5 x 10~^ respectively for t — > cZ, t ^ cry and t — > eg. 

ii) Two Higgs doublet models with and without Natural Flavor Conservation (2HDM-I, 
2HDM-II and 2HDM-III) [4], [11], [12], [13], [14], [15], [16] where the branching ratio can 
be of the order 10~^, 10"'^, 10~^ and 10^^ respectively for t cZ ^ t —> C7, t ^ eg and 

is. Note that under some assumptions on the textures of Lepton Yukawa 
couphngs, the 2HDM can also leads to Lepton Falvor violation in Higgs boson decays see 
Ref. [17] for details. 

iii) and a variety of super symmetric models [18], [19], [20], [21]. In such cases, substantial 
enhancement can be found in R parity violating MSSM models with branching ratios of 
about 10^'^, 10^^ and 10^"^ respectively for t — > cgf, t — > 07 and t eZ. The corresponding 
values in the MSSM with R parity conservation are of the order 10"^. A particular channel 
is the top FCNC coupling t ^ e^, ^ — , H° , . In the case of flavor violation induced 
by gluino, one can reach 10~^ branching ratio for t — > eh^ [20]. 
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Hence, top and Higgs FCNC ofTer a good place to search for new physics, which may 
manifest itself if those couplings are observed in future experiments such as LHC or ILC 
[1,2]. At LHC with low luminosity 10fb~^, 8 million ti pairs per experiment per year will be 
produced. This number will increase by one order of magnitude with the high luminosity 
option. Therefore, the properties of top quarks can be examined with significant precision 
at LHC. For top FCNC, it is possible to reach the following limits [1, 22]: 

Br(i cH) < 4.5 x 10"^ , Br(t ^ 07) < 3.7 x 10"^ 

Br(i ^ cZ) < 7.1 X 10"^ , Br(i ^ eg) < 10"^ (4) 

At ILC the sensitivity is less [22]:: 

Br(i ^ cH) < 4.5 x 10"^ , Br(i ^ 07) < 7.7 x 10"^ (5) 

Consequently, models which can enhance these FCNC couphngs close to the above hmits 
are welcome. 

The aim of this paper is to study top and Higgs FCNC couplings in the framework 
of the tree level Flavor Conserving two Higgs doublet Models type 1 and 11. We update 
previous results for t — >• 07, t — cZ and t ^ eg [4] and revisit the decay t — > ch^ in the 
light of full unitarity constraints and vacuum stability conditions on the Higgs sector. We 
also present our results for Higgs FCNC such as /i° — > fc, — > tc and — > tc in 2HDM 
imposing unitarity constraints. Moreover, we address the signature of those top and Higgs 
FCNC couplings at e+e^ Colliders, its 77 option as well as at future muon colliders. The 
isolated top quark signature may facilitate the search for Higgs FCNC $ — > tc events and 
also associate top-charm production. 

The paper is organized as follows. In the next section, the 2HDM is introduced. Rel- 
evant couplings are given, and theoretical and experimental constraints on 2HDM param- 
eters are discussed. In the third section, we will study the effect of 2HDM on top FCNC 
decays such as t — > 07, t — > cZ, t — > cgf and t — > ch^. Section four is devoted to Higgs FCNC 
couphngs: Br{^ ^ ic), $ = A°, which are evaluated in 2HDM-I and 2HDM-II. In 

section 5 we investigate the top-charm associate production at e+e", its 77 options and 
muon colliders. Our conclusions are given in section 6. 

2 The 2HDM 

Two Higgs Doublet Model (2HDM), is formed by adding an extra complex SU{2)l ® 
U{1)y scalar doublet to the SM Lagrangian. Motivations for such a structure include 
CP-violation in the Higgs sector [23, 24] and the fact that some models of dynamical 
electroweak symmetry breaking yields the 2HDM as their low-energy effective theory [25] . 
In particular, the Higgs sector of the Minimal Supersymmetric Standard Model (MSSM) 
takes the form of a constrained 2HDM. 
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The most general 2HDM scalar potential which is both SU{2)l ® U{1)y and CP in- 
variant is given by* [26] : 



A4(|$in$2r - \^t^2?) + A5[3?e($+$2) - v^V2f + A6[9m($+$2)]' (6) 

where $1 and $2 have weak hypercharge Y=l, vi and V2 are respectively the vacuum 
expectation values of $1 and $2 and the Aj are real-valued parameters. Note that this 
potential violates the discrete symmetry $j — > — $j only softly by the dimension two term 
A53?e($j'"$2)- The hard breaking terms (dimension four) of the discrete symmetry have 
been set to zero. 

The above scalar potential has 8 independent parameters {\i)i=i,...,&, vi and V2. After 
electroweak symmetry breaking, the combination v\ + V2 is thus fixed by the electroweak 
scale through + v| = {2\/2Gf)~^- We are left then with 7 independent parameters. 
Meanwhile, three of the eight degrees of freedom of the two Higgs doublets correspond to 
the 3 Goldstone bosons (G^, G^) and the remaining five become physical Higgs bosons: i/", 

(CP-even), (CP-odd) and H^. Their masses are obtained as usual by diagonalizing 
the mass matrix. The presence of charged Higgs bosons will give new contributions to the 
one-loop induced top and Higgs FCNC couplings (see Fig. 1 dn — > di^). 

It is possible to write the quartic coupling Aj in terms of physical scalar masses, mixing 
angles tan /3, a and A5 as follow [28] : 

9^ , ^9 , 9^ , ^9 , sin a cos a ^9 ,,9, A5 

A4 = -^Ml± , As = -^Ml , A3 = -(Ml - M?) (7) 

2m^ ^ ' ^ 2m^ ^ ' ^ Smf^. sin /3 cos ^ ^ 4 ^ ^ 

X 9'^ r 2 2 , • 2 2 sin a cos a 2 2m ^5 . , , 2 

^^ = 8cos/?^m^ ^^"^ am^ + sm am, - -^^(m^ - m,)] - -(-1 + tan /3) 

— — 2 /D — 2~N^^ am^ + cos^ CKinl — sin a cos a tan /3(m^ — m^)] + — (1 ^ 



Ssm'pmy ^ ^ '''' 4' tern' (3' 

We are free to take as 7 independent parameters (Aj)j=i 5 and tan/? or equivalcntly the 
four physical scalar masses, tan P, a and one of the A^. In what follows we will take as free 
parameters: 

A5 , MfiO , Mffo , Mao , Mh± , tan/3 and sin a (8) 

We list hereafter the Feynman rules in the general 2HDM for the trilinear scalar cou- 
plings relevant for our study. They are written in terms of the physical masses, a, jS and 
the soft breaking term A5 [28]: 

H°H+H~ = — — ^-^-—{M%o{cos (3^ sin a + sin (3^ cos a) + M%± sin 2(3 cos{(3 - a) 
M\Y sin 2(3 

— sin(/3 -I- a)X5v'^) 



*lt exist several ways of writing such a potential see [27] for a discussion 
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h^H+H" = — — (M^o (cos a cos - sin a sin + M|± sin 2/3 sin(/3 - a) 

- cos(/3 + ajAst;^) (9) 
h'H+G- = ^ cos(/3 - a){Mlo - Mf^) 



^(M^.-Ml) , .^ = ^ (10) 



We need also the couplings of the scalar bosons to a pair of fermions both in 2HDM-I and 
2HDM-II. In these couphngs, the relevant terms are as follows: 

hHto^Mt^-^ , HHto^Mt^^ , AHto,^ 2HDM-I,II 
sm p sm p tan p 

, fiT, ,,cosq; ^^nr, , ^ sin a ,nT, Mh ^ 

h%b oc Mb , H%b oc M^, , A%b oc — 2HDM - I 

sin p sm p tan p 

/i°66(xM6^^ , H%b(xMt,'^ , A°66(xMbtan/3 2HDM - II 

cos p COS p 

{H-bt)L(x-^ , {H-bt)R^-^ 2HDM-I 
tan p tan p 

{H-bt)L(x Mhtanp , {H-bt)R(x ^- 2HDM - II (11) 

tan p 

To constrain the scalar sector parameters we will use both vacuum stability conditions 
as well as tree level unitarity constraints. In our study, we use the vacuum stability 
conditions from [29, 30], which are given by: 

Ai + A3 > , A2 + A3 > , 



2^(Ai + A3)(A2 + A3) + 2A3 + A4 + min (0, A5 - A4, Ag - A4) > (12) 

It is well known that the unitarity bounds coming from a tree-level analysis put severe 
constraints on all scalar trilinear and quartic couplings [31, 32, 33]. The tree level unitarity 
bounds are derived with the help of the equivalence theorem, which itself is a high-energy 
approximation where it is assumed that the energy scale is much larger than the Z° and 
gauge-boson masses. To derive these unitarity constraints, we consider the following 
neutral scattering processes 

W^Wl, ZlZl, W^H-, WlH+, ZLh\ ZlH\ ZlA\ AM", H+R-, A'^H^ 
h%\h''H\H'^H\ (13) 

as well as the charged channels: 

Wth^, WtH^, WtA^, WtZL, ZlH+, H+h^, H+H'^, H+A^ (14) 
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It has been demonstrated, first in 2HDM with exact discrete symmetry [32] and later 
in 2HDM with softly broken discrete symmetry [33], that the above neutral scattering 
processes lead to the following 14 unitarity constraints^: 



|a±| , \b±\ , |c±| , \d±\ , |/±| , |ei,2| , |/i,2| < 167r (15) 



a± = 3(Ai + A2 + 2A3) ± V9(Ai - A2)2 + (4A3 + A4 + 0.5(A5 + Ae))^ 



with 



6± = Ai + A2 + 2A3 ± V(Ai - A2)2 + 0.25(-2A4 + A5 + Ae)^ 

c± = Ai + A2 + 2A3 ± yJ{Xi - A2)2 + 0.25(A5 - A6)2 

A5 5 A5 1 

ei = 2A3 - A4 - — + -Ae , 62 = 2A3 + A4 - — + -Ae 

5 1 11 
/+ = 2A3 — A4 + -A5 — -Ae , /- = 2A3 + A4 + -A5 — -Ae 

/i = /2 = 2A3 + ^As + ^Ae (16) 

In addition, it has been noticed in [33], that the inclusion of charged scattering processes 
leads to one extra unitarity constraint and is given by (see [33] for more details): 

pi = 2(A3 + A4) - 0.5(A5 + Ae) < IGtt (17) 

This additional constraint together with the others play an important role in constraining 
the parameter space of the 2HDM [33]. In our analysis, we take into account also the 
following constraints when the independent parameters are varied. From the theoretical 
point of view: 

i) The extra contributions to the 5p parameter from the Higgs scalars should not exceed 
the current limits from precision measurements [34]: \dp\ < 0.001. Such extra contribution 
to 6p vanishes in the limit Mh± = Ma which gives A4 = Ag. Under this constraint the 
2HDM scalar potential is 0(4) symmetric [35]. In this case {H^ , A^, H^) form a triplet 
under the residual global SU{2) of the Higgs potential. It is this residual symmetry which 
ensures that p is equal to unity at tree level. It is clear from this discussion, that 5p will 
be violated once the splitting between and is large. 

a) From the requirement of pcrturbativity for the top and bottom Yukawa couplings, tan/? 
is constrained to lie in the range 0.1 < tan/5 < 70 [36]. 

in) It has been shown in [37] that for models of the type 2HDM-II, data on B ^ Xg^ 
impose a lower limit of Mh± ^ 350 GeV. In type I 2HDM, there is no such constraint 
on the charged Higgs mass [37]. Recently, it has been shown, using a renormalization 
group improved calculation of the B branching ratio, that the scale A = m;, — £^0 

introduces significant perturbative uncertainties [38]. Such large theory uncertainties may 
have important implications for searches of new physics in radiative B decays. It is pointed 



'fWe use the rather conservative unitarity bound 167r instead of Stt 



6 




Figure 1: Generic contribution to $ — > /1/2 in SM di — > dio, in 2HDM dn dig 



out in [38] , that the lower bound on the charged Higgs mass is reduced to about 200 GeV. 
In our numerical analysis we will ignore these constraints and allow Mh± ^ 200 GeV in 
order to localize regions in parameter space where the branching ratios are sizeable. 

iv) The most recent data on the a^ = {g — 2)^/2 indicates [39] 

Aa^ = a^^P - af^ = (23.9 ± 10) x lO"'" (18) 

If the SM prediction is based on the hadronic contribution from e+e~, the above result 
(18) shows 2.4 a deviation from SM. The dominant 2HDM contribution to {g — 2)^^ come 
from one loop diagram with Higgs exchange [40, 41] and also from two-loop Barr-Zee 
diagram with heavy fermion running in the upper loop [42, 43]. In both cases, sizeable 
contribution to {g — 2)^ appear only for relatively light neutral Higgs and large tan {3. In 
our study we will assume that all neutral Higgs are heavier than 100 GeV (see below). 
This assumption make the one loop contribution to {g — 2)^ within experimental range 
for any value of tan/3 < 70 and the mixing a. The two-loop Barr-Zcc diagram has been 
studied extensively in [43] . It has been shown in [43] , that for neutral Higgs bosons heavier 
than 100 GeV and tan j3 of the order 60, {g — 2)^ two-loop Barr-Zee diagram are within 
the experimental range. Therefore, since we will assume that all Higgs are heavier than 
100 GeV we will ignore {g — 2)^ constraint in what follows. 

v) From the experimental point of view, the combined null-searches from all four CERN 
LEP collaborations derive the lower limit Mh± > 78.6 GeV (95% CL), a limit which applies 
to all models in which BR(i7^ —>■ tVt)+ BR(i7^ — > cs)=l. For the neutral Higgs bosons, 
OPAL collaboration has put a hmit on hP and masses of the 2HDM. They conclude 
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Figure 2: t — > cy, t —>■ eg and t cZ as function of tan/3 in 2HDM-I and 2HDM-II for 
Mh± = 120 GeV (left), Mh± = 250 GeV (right) 



that the regions 1 < Mh < 44 GeV and 12 < Ma ^ 56 GeV are excluded at 95% CL 
independent of a and tan j3 [44] . In what follows, we will assume that all Higgs masses are 
greater than 100 GeV. 



3 Top FCNC in 2HDM 

In this section, we will discuss the 2HDM contribution to top FCNC couplings t ^ cj, 
t ^ eg and t ^ cZ as well as t ^ ch^. In the case of t ^ cy, t ^ eg and t ^ cZ couplings, 
we reproduce and update the results of ref. [4], while for the case of t ^ ch^ couplings, we 
revisit the study done by S. Bejar et al [14] in the light of the full set of tree level unitarity 
constraints and vacuum stability conditions. 

The full one-loop calculation presented here is done in the 'tHooft gauge with the help of 
FormCalc [45] and FF packages [46]. Before presenting our numerical results, we would 
like to mention that we have adopted the following inputs: 

a = 1/137.035989 , Mw = 80.45 , Mz = 91.1875 GeV , = 1 - M^/Mj (19) 
Mt = 174.3 , Mb = 4.7 , M^ = 1.5 , M, = 0.2 GeV , = 0.04 , a,(Mt) ^ 0.105 

In the 2HDM type I and II, the top FCNC couplings t — > cy, t ^ eg and t ^ cZ are 
sensitive only to the charged Higgs mass and tan/?. We update the plots of ref. [4] for 
2 values of charged Higgs mass Mh± = 120 and 250 GeV and tan/? e [0.1,100]. Our 
results, in perfect agreement with [4], are illustrated in Fig. (2). In these plots we show 
the Branching ratios of t — > cy, t ^ eg and t ^ eZ both in 2HDM-I and 2HDM-II as 
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a function of tan/3 for Mh± = 120 GeV (left plot) and Mh± = 250 GeV (right plot). 
As can be seen from these plots, only at small tan (3 (for 2HDM-I) and large tan (for 
2HDM-II), can the branching ratios be in the range 10^*^ — 10^^. Even t eg decay which 
is proportional to the strong coupling does not exceed a branching ratio of 10^^ at large 
tan (3 > 70. Therefore, if one of those top FCNC emerge at LHC, with a rate larger than 
10~^ that would definitely attest the existence of new sources of flavor violation other than 
the one induced by the 2HDM structure. 

Let us now concentrate on i — > c/i°. This top FCNC coupling depends on the scalar 
parameters through the pure scalar couplings: h^H^H~ and h^G^H^ . These scalar cou- 
plings eqs. (10), depend on the seven scalar parameters listed above eq. (8) and may be 
sources of enhancement. Feynman diagrams which depend on these pure scalar couplings 
are depicted in Fig. (1): dyi-, dis and du. Top FCNC t — > c/i° depends also on the bottom 
Yukawa coupling through H~^bt and/or h^bb interaction which could also enhance the width 
of t — s> ch^ for large tan/9. These kinds of diagrams can be seen in di^^ and rfn^ig. 
At large tan f3, it is expected that the main enhancement comes from diagrams of the type 
du and di2. If we neglect the charm mass and right coupling [H^ ft)^ oc 1/ tan/3, which 
are suppressed for large tan /3, the amplitudes Mn of du and M12 of di2 are given by: 

Mn^^MtVftK} {H-ft)L {H-fc)L h{Ml^, M}, Ml)vt^-^u, (20) 

M12 ^ ^"felf ^i^ft^cf {H-ft)L {H-fc)L /2(Mi±,M,^M|,M,>,i^<<21) 

/i and /2 are form factors and / is the internal quark. The couplings hPH'^H~ and 

{H~ft)L are given in eqs. (10, 11). Since in 2HDM-II, behaves like m/tan/3 at large 
tan /? and the coupling h^H^H~ is enhanced for large tan (3 (see discussion below), one can 
expect that the amplitude of those diagrams grows like tan'^ f5. But care has to be taken, 
since the tan^/3 enhancement of the amplitude may be reduced by the CKM elements, 
such as Vcbi as well as by the GIM cancellation. When computing the branching ratio 
Br{t — > chP), in the top width we include both the tree level decay t — > W^h as well as 
t H^b when it is open. We make a comparison of our results with the results obtained 
by Bejar et al [14], and find perfect agreement. However the parameters set chosen by 
Ref. [14] to illustrate their results violate both unitarity and perturbativity constraints on 
Ai[32, 33]. 

We have performed a systematic scan over the full parameter in the range allowed by 
the theoretical constraints such as unitarity constraints and vacuum stability conditions 
as well as by the experimental data. As one can see from the scatter plot Fig. (4), the 
branching ratio of t — > ch^ remains less than 8 x 10~^ for light charged Higgs mass and 
large tan /3. 

Let us start by discussing the order of magnitude of the triple coupling hPH'^H~. In 
the left plot of Fig. (3) we illustrate the ratio R = \h^H+H- /h%^h\ with h%^h^ being 
the SM couphng given by h^h^h^ = 3gMlo/{2Mw)- In this plot, we use Mho = 125 GeV 
and the other parameters take the values Mj^-t = 150, M^o = 215 and Mho = 250 GeV. 
As it can be seen from this plot, unitarity constraints allow tan/3 to be as large as 70 
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Figure 3: Left, Absolute value of h^H^H coupling normalized to SM coupling h^h^h^ 
taken at M° = 125 GeV, the other parameters are fixed as Mii± = 150, M^o = 215 and 
Mho = 250 GeV, A5 = 2 . Right, 2 x Br{t ch°) in 2HDM-1I as function of tan/? 
for several values of Mh± and M^o = 125, M^o = 215, M^o = 250 GeV, A5 = 2 and 
sin a = 0.29 
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Figure 4: Scatter plots for Brit ch!^) in the (tan/?, Mh±) plan (left) and (M/jO, Mj 
plan. We scan over: 100 < M?,o < 175, 100 < Mho,ao,h± ^ 600 GeV, 
20 < tan/5 < 70 and A5 = Ml/v"^. 



-1 < sin a < 1 and 



only for sin a ^ 0.29 and A5 = 2. For large sin a = 0.8 unitarity already breaks down 
for tan/? ?s 10. It is clear that the coupling h^H^H~ reach its maximal value for large 
tan /3. The minimum observed in the h^H^H~ coupling is due to a cancellation between 
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the different terms which contribute to hPH^H^ (eq.(9)). Such a minimum depends of 
course on the values of A5 and sin a. Since the couphng h^H^H~ depends hnearly on A5, 

it is obvious that this couphng gets its largest value for the largest value of A5 tolerated by 
perturbativity: A5 ~ Stt. However, unitarity requires that for large A5, large tan/? is not 
allowed (see left plot Fig. 3 for A5 = 6). 

Having the order of magnitude of hPH'^H~ in mind, we show in the right panel of 
Fig. (3) the Branching ratio ol t ^ chP as a function of tan /3 for several values of Mh± 
with Mfeo = 125, Mao = 215, M^^o = 250 GeV, A5 = 2 and sin a = 0.29. For large tan p and 
light charged Higgs Mh± = 150 GeV, Br(t ch^) can reach value larger than 10~^ while 
for Mh± = 200 GeV, the Br(t ch^) is already less than 10~^ even for large tan/9 ~ 70. 
This enhancement is mainly attributed to the triple scalar coupling h^H^H^ as well as to 
the bottom Yukawa coupling which are both enhanced for large tan (5 as one can see from 
eqs.(20,21). 

We also provide scatter plots Fig. (4) for Br(t — ch^) in {Mh±, tan (3) and {Mh±, M^o) 
plans. There is only a small window for large tan j3 and Mh± < 180 GeV where Br(t ch^) 
can be larger than 10~^. As argued before, due to the large theoretical uncertainties in 
the B Xgj calculation, the lower bound on the charged Higgs mass is reduced to about 
200 GeV [38]. From the left panel of Fig. 4, it is illustrated that for charged Higgs mass 
Mh± e [200, 300] GeV and large tan/?, Br(t ch^) G [10^^ 10"^]. One concludes that the 
top FCNC t ch^ in 2HDM can reach a detectable rate at LHC with the high luminosity 
option. 

4 Higgs FCNC in 2HDM 
4.1 Higgs FCNC in SM 

Before presenting our results in 2HDM, we would like to comment on the Branching ra- 
tio of if — > fc in the SM. Recently, an estimation using dimensional analysis and power 
counting, appeared both for Br{H — > sb) [47] and Br{H — > tc) [15] while in [16], the full 
diagrammatic calculation has been done. Here we comment on results for Br{H — > tc) and 

Br{H 6s)[16]. 

As expected the branching ratios of these decays are very suppressed due to GIM mecha- 
nism [15, 16, 47]. In the case of if — > tc, the internal fermions are bottom, strange and down 
quarks with masses close to degeneracy. The GIM suppression in this case is very severe, 
and the Br is of the order 10"''^^. For if — > 6s the internal fermions are top, charm and up 
quarks and with the large top mass as well as large splitting between the internal quarks, 
the Branching ratio of if 6s is in the range 10~^-10~^. if — > 6s has another advantage 
over if — > fc which is the fact that if — > 6s can be open for Higgs mass Mh < 2Mw and 
hence the decay width of the Higgs can be very narrow which can enhance the Branching 
ratio for Mh ^ 2Mw- The decay if ° — > ic is open only when Mh > rrit + rric ~ 2Mz- For 
a Higgs mass in this range, the decay channels if° — > {W~^W' , Z^Z^} are already open 
and so the total width of the Higgs is no longer narrow. Such a large decay width may 
reduce the Branching ratio of if° — > tc. For instance, for Higgs mass heavier than 250 
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GeV, we get a branching ratio of the order 10 — > 10 for H — > tc. New sources of 
flavor violation hke charged Higgs are expected to enhance the branching ratio of $ — > fc. 



4.2 $ ^ 

In this section we discuss the 2HDM contributions to the Higgs FCNC. We concentrate 
mainly on hP tc, — > tc and — > tc. The generic Feynman diagrams which contribute 
to those FCNC are shown in fig .(1). As wc can sec. the additional interactions like 
^^H+H^, ^^H+G~ and give new contributions to $° ic. 

The calculation of — > tc, — > tc and to in 2HDM has been first done in 

Ref. [15]. In Ref. [15], only trilinear Higgs self-couphngs Xhhh have been constrained by 
the maximum unitarity limit tolerated for the SM trilinear coupling: 

lA^^^I < iXfTHHiMH = lTeV)| = ^^^1^ (22) 

We have checked that the above constraints (22) are not enough to guarantee unitarity 
and perturbativity constraints. For instance, for the parameter set considered in Ref. [15], 
Ai can take extremely large values of the order ^ 5000 which is not tolerated neither by 
perturbativity nor by unitarity constraints. We perform a cross check with Ref. [15], and 
our results^ are about 3 times larger than those presented in Ref. [15]. 

The decay width r$ of all scalar particles h^, and AP are computed at tree level in 
2HDM as follows: 

r$ = ^ r($ ^ //) + r($ ^ vv) + r($ ^ VHi) + r($ ^ HiHj) (23) 
/ 

QCD corrections to $ ^ // and $ {gg,j^,jZ,V*V*,VV*,V*Hi} decays are not 
included in the width. Since the width of $ ^ tc does not become comparable to $ — > 
{gg, 77, '~fZ}, it is not necessary to include them in the computation. The full width of the 
Higgs bosons is taken from [48] without QCD corrections. As stated before, we perform a 
systematic scan over the 2HDM parameter space using unitarity constraints and vacuum 
stability conditions on the scalar sector parameters. In our numerical analysis, we limit 
ourself to 2HDM type II. As shown in Ref. [15], the rates in 2HDM-I are smaller than 
the rates in 2HDM-II. We also present numerics for CP even Higgs only $ = and 
comment on our finding for — > tc in 2HDM type I and II. 



h^^tc 

We illustrate first in Fig. (5) the branching ratio of — > tc in 2HDM-II as a function of 
Mj^o (left) and tan/3 (right). The parameters are fixed as: M40 = 215, Mjjo — 250 GeV, 
As = 2 and sin a = 0.5 and Mh± = 150 GeV (left) and M^o = 220 GeV (right). In both 
plots (left and right), one can see that the Branching ratio of — > tc, like the Br(t — > ch^), 

■tWe reproduce our results by 2 different and independent programs. 
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Figure 5: 2 x Br{h^ — >• tc) in 2HDM-II as a function of Mh (left) and tan (3 (right) with 
Mao = 215, Mho = 250 GeV, Ag = 2 and sin a = 0.5 and Mh± = 150 GeV (left) and 
MhO = 220 GeV(right) 




Mh+ (GeV) Mh+ (GeV) 

Figure 6: Scatter plots for 2 x Br(/i° tc) in the (tan/3, Mh±) (left) and (M^o, Mh±) 
plans. We scan over: 180 < M/,o < 230, 190 < Mho < 230 100 < Mao,h± ^ 600 GeV, 
-1 < sina < 1 and 20 < tan/5 < 70 and A5 = M^/v'^. 

can be of the order 10~^ only for large tan/5 and relatively light charged Higgs mass. As 
in the case of t ^ ch^, this enhancement at large tan/3 is due to bottom Yukawa coupling 
as well as to the triple scalar coupling h^H^H^ . 

We show on the right panel of Fig. (5) the Branching ratio of — > tc as a function of 
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Figure 7: Branching ratio for 2 x Br{H^ tc) as a function of Mho (left) and tan/5 (right) 
in 2HDM-II Mho = 120, and Mao = 200 GeV, A5 = 1, sin a = 0.75 with Mh± = 110 GeV 
for left and M^o = 220 GeV for right panel 



tan/3 for several values of charged Higgs mass Mh±- The maximum rate for Br(/i° 
tc) ~ 4 X 10"^ is obtained for Mh± = 110 GeV and large tan/? 65. The reason for this 
enhancement is due to threshold effect, since at M/^o = 220 GeV and Mh± = 110 GeV the 
decay channel H^H~ is open^. For charged Higgs mass away from the threshold 

region, the 'Bi{hP — > ic) is reduced and for Mh± = 150 GeV, it is of the order 10~^ at 
tan 13 ^ 70. 

In the scatter plots Fig. (6) it is illustrated in the left panel that for large tan/3 > 30 
and light charged Higgs Mh± < 250 GeV the Br(/i° —>■ ic) G [10~^, 10~^]. The maximum 
branching ratio is obtained close to threshold region where the decay —>■ H^H~ is open. 
For charged Higgs mass Mh± ^ 250 GeV, Br(/i'^ ic) < 10~^ for all range of tan/?. From 
the right panel, one can see that the Br(/i° — > ic) is not very sensitive to the CP-even mass 
Mho. 



^tc 

We now discuss the 2HDM contribution to —* ic. We perform a scan in the 2HDM 
parameter space and find a region where the branching ratio is sizeable. The conclusion 
is similar to —>■ ic. To have a large branching ratio we need to maximize the coupling 
H^H^H~ and the bottom Yukawa coupling. This is achieved by taking large tan (3. 
In the left panel of Fig. (7), we illustrate the Branching ratio of — > tc as a function of CP 
even Higgs mass Mh for several values of tan (3. As it can be seen, for moderate tan /3 < 15, 
the Branching ratio is less than 10~^. When the CP even mass Mho ~ 2Mh± ~ 220 GeV 

^see next plot for similar threshold effect for H'^ ic 
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Figure 8: Scatter plots for 2 x Bi{H° ic) in the (tan/3, Mh±) (left) and (M^^o, Mh±) 
plans. We scan over: 100 < M/,o < 230, 180 < Mj^o < 600 100 < Ma'-\h± ^ 600 GeV, 
-1 < sina < 1 and 20 < tan/3 < 70 and A5 = Ml/v"^. 



one can see a peak. Such a peak corresponds to the opening of the decay channel 
H^H^ . This threshold effect is also amplified by the fact that the coupling H^H^H^ 
takes its maximal value for Mho ~ 220 GeV and large tan (3. In this region, the Branching 
ratio can reach 3 x 10~^ at large tan/3 50. As can be seen in the left plot, for large 
tan/3 > 25 unitarity is violated apart from a small region around Mh ~ 220 GeV. Once 
the Higgs mass Mh becomes larger than 2Mh±i the decays —>■ {H^H~ , h^h^} are open, 
the width Fj^o becomes large and hence the Branching ratio is reduced. 

In the right panel of Fig. (7), we show the branching ratio for —>■ tc as a. function of 
tan/3 for several values of Mfj±- As for the — > tc case, the branching ratio is enhanced 
by large tan/3 and rather light charged Higgs mass. Also in this case, for large tan/3 
the coupling H^H^H~ reaches its maximal value. It is clear from the right plot that 
the Branching ratio Br{H^ tc) can be greater than 10~^ only for large tan/3 > 45 and 
Mh± = 150 GeV. Once the charged Higgs mass becomes large, the Br{H^ tc) decreases. 
This can be seen in the right plot of Fig. (7) where Mh± = 250 GeV, the Br(if° — > tc) is 
less than 10~^. 

We have also performed a systematic scan for this case, and the results are illustrated in 
Fig. 8. In those scatter plots, we show Br(if° — > tc) as function of (Mff±,tan/3) (left) and 
[Mh±,Mho) (right). Again, close to the threshold region where —>■ H^H' is open, 
the Br(if° ic) can be larger than 10~^ (see left panel of Fig. 8). Such an enhancement 
requires of course Mh± ^ 150 GeV and large tan/3. In this case also, there is a region in 
which the charged Higgs mass is larger than 200 GeV, which could accommodate B — > ^^7, 
and the Br(i^'^ — > ic) G [10~^, 10~^]. From the right panel of Fig. 8, one can see that once 
the it threshold has been passed Mh ^ 350 GeV, the Branching ratio Br(if° ic) is 
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Figure 9: Topological contributions to e^e — > tc and jji^jj, — > tc in 2HDM 



reduced. 

^tc 

At the end of this section, we would like to comment on tc decay in 2HDM type 1 and 

II. In 2HDM-II, we perform a systematic scan over the 2HDM parameters. It is found that 
the Branching ratio of — > tc decay never exceeds 10^* (resp 3 x 10~^) for small (resp 
large) tan/3. The reason is that, unlike hP and the CP-odd ^4°, being pseudo-scalar, 
does not couple to a pair of charged Higgs. Therefore A'^ — > tc, has only a unique source 
of enhancement which is the bottom Yukawa coupling. 

In 2HDM-I, there is a smaU region for 100 < Mh± < 200 GeV, 176 < M^o < 2Mt ^ 350 
GeV and small 0.1 < tan/3 < 0.3 where the Branching ratio of A^ — > tc is in the range: 
[10~^,3 X 10^^]. At large tan/3 the branching ratio decreases quickly and is less than 
10~^^ for tan/? > 2. The origin of this is the fact that the couplings A^bb and H~tb are 
proportional to 1/ tan/3 and so suppressed for large tan/3. 



5 Top-charm production at e+e , 77 and muon collid- 
ers 

There have been several studies looking for collider signatures of the top and Higgs FCNC 
couplings both at lepton colliders [49] as well as at hadron colliders [50] . Because of the 
very large mass of the top. Top-charm production at leptons and hadron colliders would 
have a clear signature. 

At LHC, the production and subsequent decay of the neutral CP-even Higgs mediated by 
FCNC have been analyzed in [15] for 2HDM and [50] for MSSM. It has been demonstrated 
in those studies that the production rates for the bs + bs channel can be large while the 
rates for tc -\- tc channel are more modest. Detection for the former channel at LHC is 
difficult but it is easier for the later [50] . 

Associate Top-charm production at e+e" colliders has been studied in SM [51] and found to 
be very small. New physics contribution to e+e^ tc such as SUSY can enhance the cross 
sections by several orders of magnitude with respect to SM [52] . The top-charm production 
at lepton colliders e'^e" and has been studied in two Higgs doublet models without 
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Natural Flavor Conservation 2HDM-III which can lead to measurable effects [49, 53]. It 
has been pointed out that the tree level vertex $tc can be better probed through t-channel 
W-W and/or Z-Z fusion at high energy e+e" collisions e+e^ — > tcviVe and e^e~ — > tce^e" 
[54, 55]. An interesting feature of these reactions is that, being t-channel fusion processes, 
their cross sections grow with energy unlike s-channel reactions e^e~ — > tc which are 
suppressed at high energies. The cross sections of e+e~ — ^ tcvj>^ and e"*"e~ tce~^e~ are 
found to be one or two orders of magnitude bigger than the cross sections of e+e" tc 
[54, 55]. 

In the present study, we limit ourself to e+e" colliders, its 77 option and muon colliders. 
The cross sections for e+e~ — > tc, 77 — > tc and/or h'^/jl" — > tc in 2HDM, if sizeable, can 
give information on the top FCNC couplings t —>■ cj and t — > cZ as well as Higgs FCNC 
$ to. Unlike the process e~^e~ to, which is s-channel suppressed at high energy, the 
process 77 —>■ tc has a t and u-channel contributions, its cross section may grow with the 
center of mass energy. 

At e+e~ coUisions, the generic contributions to e'^e~ — > tc are drawn in Fig. (9). We 
have three sets of diagrams: 'jtc vertex contributions through e^e~ — j> 7* — ;> to, Ztc 
vertex contributions through e~^e~ Z* ^ tc and box contributions. The topological 
contributions to e^e~ — > ic are depicted in Fig .(9). 

In 77 collisions, we have more topological diagrams that contribute to 77 ic. jic vertex 
enters game through diagrams d2, ds, and ds of Fig. (10) and their crossed one. The 
quartic coupling of 2 photons to a pair of H^H~ or W^W~ gives the topological diagram 
di. Box (resp self-energies) topologies are depicted in diagrams ^4 and d^ (resp d^). We 
have checked both analytically and numerically that the result is Ultra- Violet finite and 
renormalization scale independent. In the case of 77 — > tc, gauge invariance has been 
checked numerically. 

Due to GIM mechanism, in the SM, the cross sections of e+e^ ^ ic and 77 ic are 
very suppressed see Fig. (11) (left and right). It is clear from Fig. (11) that, in the SM, 
the 77 ic cross section is about one order of magnitude bigger than e+e" ic cross 
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Figure 11: Cross sections for associated top-charm production at e~^e~ colliders (left) and 
its 77 option (right) in 2HDM-II with tan/3 = 65 and several choices of charged Higgs 
mass 



sections. Our results for the SM agree with [51] 

We illustrate our numerical results in 2HDM-II for large tan (3 = Q5 and several values of 
charged Higgs. The reason for this choice is the fact that the one loop effective coupling 
t ^ cj and t ^ cZ are larger for large tan/3 (see Fig. (2)). We consider only unpolarized 
cross sections. It is well known that the cross sections for polarized initial states differ from 
the unpolarized cross sections only by the normalization factor. 

In both cases, for e~^e~ and 77 collisions, the cross sections in the 2HDM are enhanced 
by about four orders of magnitude with respect to their SM values. Above y/s = 200 GeV, 
the cross sections become sizeable and develop a peak at ^/s = 2Mh±, which is due to the 
threshold effect of charged Higgs pair production. The cross sections can reach a value of 
0.003 fb (resp 0.1 fb) for e+e~ — > ic (resp 77 — > ic) near threshold. At 77 options, cross 
sections of about 0.01 — 0.1 fb would give 4 — > 80 events for the assumed luminosities of 
400 — > 800 fb~^. Above the threshold region i/i > 800 GeV, the cross section decreases 
with increasing the energy and reaches a value of ~ 10~^ fb at ^/s > 1.5 TeV. 
In the 2HDM-I, the cross sections are below ^ 10^^ fb (resp 10~^ fb) for 77 (resp e^e~) 
collisions. 

In recent years an increasing amount of work has been dedicated to the physics pos- 
sibilities of /i"*"/!" colliders [3]. Such colliders offer novel ways of producing Higgs bosons, 
such as an s-channel resonance in the case of neutral scalars. The production process 
/i+yU" ic proceed through the topological diagrams drawn in Fig. 9. As for e+e^ ic 
we have the gauge boson s-channel ^~ 7, Z* — > ic contributions as well as the Higgs 
boson s-channel exchange [i~ —>■ h*, H*, A* — >• ic. The three couplings hP^i^ ii~ , H'^fi~^fi~ 
and A^fi~^fi~ behave like tan /3 at large tan f5 and this might be a source of enhancement 
in the large tan/3 limit. Such enhancement has been already observed for /x^/x^ W^H^ 
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Figure 12: tc + tc cross section as function of ^/s (left) and tan/? (right) for 

Mh = 100,Mh = 278, Ma = 170 GeV, sin a = -0.84 and A5 = M^/v^ = 0.95. 

and ^ ZA^ [56]. 

In the SM, it is expected that the cross section /i"*"/!" — > ic is very suppressed as in 
the case for e"'"e~ tc. In the case of fi^fi~ tc, the only difference with respect to 
e~^e~ — > ic is the contribution of Higgs boson s-channel H* tc. However, in SM 

the couplings H^fi~^fi~ as well as H^ic does not have any enhancement factor like tan/3. 
Therefore, the cross section for fi^^^ ic will be similar in magnitude to the e"'"e~ ic 
cross section. 

We show in Fig. (12) cross sections for 11^ ^~ ^ tc as a function of center of mass energy 
(left) and as a function of tan/3 (right). As expected, the cross sections are enhanced for 
large tan/3 where both 74°}/i+/i~ and A°}tc couplings are enhanced. As 

can be seen from right plot of Fig. 12, the cross section can reach a few fb near threshold 
region for tan /3 = 65 and light charged Higgs. One can see also several kinks, the first one 
on the left is located around y/s ~ Mh = 278 GeV and is due to s-channel resonance of the 
heavy CP-even Higgs . The second kink, located at several points where ^/s ~ 2Mh±, 
is due to the threshold production of pair of charged Higgs. For Fig. 12, instead of using 
unitarity constraints we have used the perturbativity constraints on the Xf Aj < 87r. That 
has been done in order to allow a larger reach for tan (3. For instance, for charged Higgs 
of 150 GeV, unitarity constraints are satisfied only for tan/3 < 48 while perturbativity 
constraints on the Aj are satisfied even for tan/3 ~ 100. As it can be seen from those plots, 
the cross section can be larger than 0.01 fb only for low center of mass energy ^/s < 500 
GeV and rather light charged Higgs. For high center of mass energy ^/s > 500 GeV, the 
cross sections of the s-channel process ic are reduced. 

We note that in 2HDM-I, ^~ — > ic would not register any observable signal. This is due 
to the fact that some fermionic couplings are proportional to cot^ /3, and so unacceptably 
small values of tan /? would be required in order to allow observable cross-sections. 
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6 Conclusions 



In the framework of 2HDM with natural flavor conservation, we have studied various top 
and Higgs flavor changing neutral couplings: t 07, cZ, eg, ch^ and $ ic. The study 
has been carried out taking into account the experimental constraint on the p parameter, 
unitarity constraints, and vacuum stability conditions on all scalar quartic couplings Aj. 
Numerical results for the branching ratios have been presented. We emphasize the effect 
coming from bottom Yukawa couplings and pure trilinear scalar couplings such as hPH^H~ 
and H^H+H-. 

We have shown that, in 2HDM-II, the branching ratios of top FCNC decays t — > cj, cZ, eg, chP 
as well as Higgs FCNC {hP, H^} tc can be orders of magnitude larger than their SM 
values. Such enhancement requires large tan/3, rather light charged Higgs and not too 
heavy CP-even Higgs masses. We have shown that t ch^ as well as {h^ , H^} tc 
can reach observable rates at LHC with Branching ratios of the order 10~^ while all the 
unitarity and vacuum stability constraints are fulfilled. With Branching ratio of this size, 
the production and subsequent decay of the neutral CP-even Higgs mediated by FCNC 
can lead to a few hundred events at LHC with the High luminosity option [15] . 
We have also studied the top-charm associated production at e+e~ colliders, its 77 option 
and also at muon colliders. It has been shown that, in 2HDM-II, the cross section of 
77 — > tc can be of the order 0.01 — > 0.1 fb near the threshold region while the correspond- 
ing cross sections at e'^e" remain well below 10~^ fb. At muon colliders, the situation is 
slightly better due to Higgs boson and A° s-channel contribution. We showed that 

//"•"/i" —>■ tc production offers an attractive way of searching for FCNC at such colliders. 
The cross section grows with increasing tan/3 with values of few fb being attainable for 
tan (3 > 48. 
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